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ALUMINIZED  ORGANIC  EXPLOSIVES 


Prepared  by: 
Donna  Price 


ABSTRACT:  Available  detonation  velocity  (D)  and  pressure  data  for 
aluminized  organic  explosives  have  been  reviewed.  The  observed 
decrease  in  D  caused  by  addition  of  aluminum  can  be  explained  theo¬ 
retically  on  the  basis  of  entrainment  and  a  small  amount  of 
compression  of  the  Al.  But  the  large  particle  size  effect  and  tne 
porosity  effect  cannot  be  so  explained.  It  has  been  suggested  that 
an  endothermic  reaction  could  explain  the  enhanced  effect  on  D  with 
decreased  particle  size  of  the  Al.  Consideration  of  recent  studies, 
however,  leads  to  the  conclusion  thav  Al  is  chemically  inert  during 
detonation.  Hence,  the  particle  size  effect  and  also  the  porosity 
effect  seem  better  explained  by  heat  transfer  to  the  Al. 
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|  This  work  was  carried  out  under  tasks  MAT  03L  000/R011  01  01  FR  59 
|  and  ORDTASK  033  102  F009  06  01  which  has  now  been  terminated.  It 
1  was  done  to  provide  a  comparison  between  aluminized  organic  explo¬ 
sives  and  aluminized  ammonium  perchlorate,  the  subject  of  a  special 
study  described  in  a  parallel  report.  The  present  results  and 
conclusions  on  aluminized  H.E.  should  be  of  general  interest  in  the 
area  of  explosive  effects. 
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INTRODUCTION 


About  13  years  ago  the  writer  reviewed  the  literature  and 
available  data  for  aluminized  organic  H.E.  At  that  time  the  various 
investigators  agreed  on  one  point,  namely,  that  aluminum  lowers  the 
detonation  velocity  (D)  of  H.E.  On  the  question  of  how  it  caused 
the  lowering,  there  was  violent  disagreement.  It  is  the  purpose  of 
the  present  report  to  review  more  recent  data  and  to  reassess  our 
knowledge  of  aluminized  explosives  with  particular  emphasis  on  the 
effect  on  D  of  adding  A1  to  organic  H.E. 


A  most  helpful  current  publication1  presents  data  on  the  effect 
of  adding  0  to  15%  by  volume  of  5u  A1  to  HMX/Viton,  85/15.  The 
motion  of  the  metal  walls  in  the  cylinder  test  was  used  to  compute 
the  energy  release  of  the  test  explosive  relative  to  that  of  HMX. 

The  most  important  contribution  of  reference  (1)  to  this  review  is 
that  the  reaction  of  aluminum  with  the  detonation  products  did  con¬ 
tribute  to  the  energy  release  but  not  until  4  usee  after  the 
detonation.  Its  maximum  contribution  was  reached  at  about  13  usee  at 
an  energy  indicating  complete  reaction  (by  relative  Ruby  cod^  results) 
of  half  the  A1  present  in  a  one-inch  dia  cylinder  test.  Thi  .  seems 
strong  evidence  that  finely  divided  A1  does  not  react  chemically  with 
the  detonation  products  of  organic  explosives  in  time  to  affect  D, 
i.e.,  in  times  of  0.1  usee  and  less.  It  is  not  conclusive  be  use 
there  is  the  possibility  of  a  rapid  endothermic  reaction.  However, 
equilibrium  computations  have  not,  as  yet,  revealed  such  an  endo¬ 
thermic  reaction  that  is  thermodynamically  possible. 


It  was,  and  is,  generally  agreed  that  admixed  aluminum  in 
sufficiently  fine  form  will  react  exothermally  in  time  to  contribute 
to  some  effects.  Inasmuch  as  the  time  available  for  contribution  in 
the  NOL  fragmentation  test  is  about  15  usee2  ,  5u  A1  should  have  a 
significant  effect  on  the  fragmenting  ability  of  charges  in  that 
particular  test  configuration.  For  blast  and  underwater  damage, 
the  available  time  for  reaction  is  one  to  two  orders  of  magnitude 
greater  for  one-pound  charges,  and  even  larger  for  larger  charges. 

It  is  common  knowledge  that  even  coarse  A1  can  increase  the  damage 
in  such  cases. 


EFFECT  OF  ALUMINUM  ON  DETONATION  VELOCITY 


On  the  evidence  given  above,  it  seems  likely  that  the  A1  in  an 
organic  H.E.  is  an  inert  diluent  as  far  as  its  effect  on  D  is 
concerned.  Since  it  does  lower  D,  it  must  do  so  by  decreasing  the 
energy  available  for  propagating  the  front  by  one  or  more  of  the 
following  mechanisms: 
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a.  entrainment  in  the  gas  products, 

b.  absorption  of  heat, 


c.  compression. 

In  an  earlier  study  of  solid  diluents  in  H.E.3,  we  found  that  air 
(or  voids)  was  the  most  efficient  diluent  available  for  lowering  D. 

To  take  care  of  this  very  large  effect  in  porous  charges,  valid 
comparison  of  solid  diluents  must  be  made  at  equal  volume  fractions 
of  air  (ya).  Comparisons  at  the  same  absolute  density,  although  still 
popular,  are  invalid. 

Thus,  examination  of  NOL  data  for  aluminized  TNT,  RDX ,  and 
TNETB4  revealed  that  the  void  content  of  the  charge  was  of  comparable 
importance  with  the  A1  content  in  determining  the  extent  to  which  D 
is  lowered.  Hence,  all  curves  here  are  given  in  the  form 


D(A,u)  =  a(u)  +  b(u)A  (] 

where  A  is  the  ratio  PQ/PV  (actual  density  to  voidless  density)  and 
u  =  d-1  is  the  reciprocal  diameter.  Comparisons  are  then  carried 
out  at  the  same  A  value . 


The  best  measure  of  the  effect  of  A1  on  D  would  be  given  by  the 
infinite  diameter  results: 

Di(A)  =  aQ  +  b0A  (2) 

where  aQ  and  bQ  are  constants  and  specific  to  the  chemistry  of  the 
explosive  tested.  Unfortunately,  the  few  diameter  studies*  available 
produce  data  far  toe  sparse  to  measure  the  effect  of  Al.  Conse¬ 
quently  ws  have  used  the  NOL  data4,  obtained  at  a  constant  dia  of 
50.8  nun,  and  have  rejected  all  data  for  solid  diluent  contents  above 
about  32  weight  %.  There  are  two  reasons  for  this.  First,  any 
diameter  effect  on  the  measured  rate  would  be  expected  to  increase 
with  increasing  amounts  of  an  inert  diluent.  Second,  several  of  the 
curves,  b(u)  vs  %  Al ,  show  a  sharp  change  of  slope  at  this  point.4 
It  is  hoped  that  the  trends  obtained  at  d  =  50.8  mm  will  not  be  too 
distorted  by  the  diameter  effect  and  hence  will  not  differ  greatly 
from  those  that  would  be  obtained  at  u  =  0 ,  provided  %  diluent  £  32%. 

TNT  Mixtures 


*0ne  on  HBX-1  and  one  on  a  tri tonal. 


Figure  la  shows  the  effect  of  Al  on  lowering  the  detonation 
velocity  of  TNT  as  D  vs  A  curves  for  0  to  32%  Al .  In  this  case, 
all  charges  had  a  diameter  of  50.8  mm  and  the  Al  was  Navy  Spec, 
atomized  material  which  passed  through  a  No.  100  sieve  (149y  opening).4 
Figure  lb  shows  analogous  curves  for  the  effect  of  NaCl  as  a  diluent. 
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In  these  charges,  the  NaCl  was  -20  +  28  mesh  (840  -  590p)  and  the 
charge  diameter  was  unspecified  although  it  was  intended  to  be 
large  enough  to  produce  infinite  diameter  D  values.5  Figure  1 
compares  the  effect  of  the  two  diluents  at  the  same  A  and  at  various 
weight  fractions  of  the  d.  a;ent.  A  somewhat  better  comparison  is 
on  an  equal  volume  basis* .  this  assures  that  diluents  will  be  compared 
at  the  same  ye  (vol  me  fraction)  of  the  energy  producing  organic 
explosive.  The  diluent  effect  on  D  can  be  expressed  by  the  ratio 
D/Da  where  Da  is  the  velocity  of  the  undiluted  TNT  at  the  same  A 
value.  Figure  2  shows  curve.,  of  D/Da  vs  yd  (volume  fraction  of 
diluent)  at  A  values  of  0.95  and  0.80.  It  is  quite  evident  that  at 
equal  volun.-  fractions,  salt  lowers  ".‘he  rate  more  than  A1  if  the 
charge  is  r  ar  vcidless  density,  whereas  the  roles  of  the  diluents 
are  revers  for  .nor*  porous  charges.  Neglect  of  this  large  porosity 
'S  much  of  the  conf  ision  in  the  older  literature  on 
which  a-.  :i+‘.on  of  A1  affects  D.  Both  Figures  1  and  2 
li  1  -nt  =:.i  t,  <.th'  lowering  of  L)  increases  as  the 
lases;  .  '  ,  therefore,  easier  to  measure  at  ?ow  A. 


effect  acc 
the  extent  .• 
show  that  th 
porosity  inc 


In  addition  to  tne  compaction  and  the  diameter  effect,  one  other 
factor  that  should  be  controlled  for  va? id  results  is  the  particle 
size  and  shape.  The  effect  of  this  last  factor  was  not  studied  in 
the  initial  NOL  jrk.  Its  importance  is  demonstrated,  however,  in 
data  published  by  Dremin  et  al.6  These  workers  were  well  aware  of 
the  importance  of  po?  >sity,  and  made  their  measurements  at  85.0  and 
96.5%  TMD.  The  D  vs  A  lines,  each  constructed  from  the  reported  pairs 
of  experimental  points,  are  shown  in  Figure  3  in  which  the  TNT/A1, 
85/15,  curves  are  very  similar  to  those  of  Figure  la.  But  Figure  3 
supplies  the  additional  information  that  ?70u  Al  has  little  or  no 
effect  on  D  whereas  0.2p  Al  has  a  very  large  effect.  (There  was  also 
one  measurement  made  with  8 Op  Al  at  85%  TMD;  that  showed  an  inter¬ 
mediate  effect.) 


If  the  NOL  data  of  Figure  la  are  interpolated  to  TNT/A1,  85/15, 
they  give  a  value  of  6.10  mm/psec  at  A  =  0.85.  Moreover,  the  data  of 
reference  (5),  when  plotted  D  vs  6(average  particle  size)  or,  better, 
D  vs  log  6,  as  in  Figure  4,  give  reasonably  smooth  curves.  Entering 
them  at  the  NOL  value  of  D  for  this  composition  (crosses  on  Figure  4) 
shows  corresponding  6  values  of  22  and  18p  respectively.  On  this 
basis,  the  Al  used  in  the  NOL  work  appears  to  have  had  an  effective 
diameter  of  about  20p  on  the  scale  of  reference  (6).  In  this  case, 
the  charge  diameter  is  not  given;  hence,  the  particle  size  scale  may 
be  influenced  by  a  charge  diameter  effect. 


To  return  to  Figure  3,  the  curve  for  85/15,  TNT/fine  W,  Is 
also  shown  from  the  reference  (6)  data.  In  this  case,  extrapolation 
to  A  =  1  does  not  give  the  100/0  value  of  D,  as  seems  to  be  the 
case  for  TNT/A1,  but  falls  somewhat  below  it  (by  less  than  3%).  It 
is  similar  to  results  of  the  same  extrapolation  for  NaCl;  in  Figure  lb 
the  curve  for  TNT/NaCl,  70/30,  extrapolates  to  about  6%  below  the  TNT 
value  at  A  =  1. 
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FIG.  3  EFFECT  OF  15%  DILUENT  ON  DETONATION  VELOCITY  OF  TNT 
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Dremin  et  al .  also  reported  results  from  charges  containing  15%  of 
fine  and  of  coarse  SiC>2  •  At  85%  TMD,  these  lie  above  the  TNT/G.  2y 
Al  curve;  at  96.5%  TMD,  below  it.  The  sharp  change  in  the  effect 
of  this  diluent  was  attributed  to  a  phase  change  ir.  the  Si02  .*  Hence* 
this  component  is  not  comparable  to  the  others  (Al,  W,  NaCI)  which 
show  no  sucn  anomalous  effect.  Certainly  the  data  cannot  be  used 
to  construct  a  D-A  curve  if  the  diluent  is  in  different  phases  at  the 
two  A  values.  Consequently ,  no  curve  for  TNT/Si02  ,  85/15,  is  given  in 
Figure  4. 


Reference  (3)  reviewed  the  attempts  to  predict  D  of  organic 
explosives  to  which  inert  solid  diluents  had  been  added.  None  of 
the  models  was  completely  successful  and  only  one  (that  considering 
heat  transfer)  incorporated  the  effect  of  particle  size.  The  work 
did,  however,  give  guidance  in  the  dependence  of  D/D„  on  volume 
fraction  of  the  diluent  (y^),  density  of  diluenx  (Pd)  and  compressi¬ 
bility  of  diluent  (a).  Patterson7  assumed  the  diluent  chemically 
inert  but  entrained  in  the  detonation  reaction  products.  He  then 
derived  the  relation 


£  S 


g 


D/Da  =  /P77y 


yepv 


(3) 


where  ye  is  the  volume  fraction  of  the  explosive,  Pe  the  density  of 
the  explosive  (1.654g/cc  for  TNT)  and  Pv  the  density  of  the  voidless 
mixture,  H.E. 'diluent .  On  the  other  hand,  Copp  6  Ubbelohde*  assumed 
that  the  inert  diluent  is  not  entrained  but  is  heated  both  by  com¬ 
pression  and  heat  conduction  from  the  hot  detonation  products.  They 
felt  that  their  experimental  data  (D  measured  on  heavily  cased 
charges  of  TNT  mixtures  including  TNT/A1)  were  satisfactorily 
reproduced  by  their  computations  based  on  these  assumptions.  To 
compute  D,  they  used  Marshall's  equation 


D(mm/ysec)  =  0.430  /nT  +  3.5(P0-1) 


(4) 


where  n  is  the  number  of  moles  of  gas  produced  per  gram  of  explosive, 
PQ  is  the  loading  density,  and  T  is  the  effective  detonation 
temperature.  Obviously  a  good  many  other  assumptions  had  to  be 
made  to  compute  T;  the  final  results  gave  the  right  order  of  magni¬ 
tude  for  the  diluent  compressibilities,  but  sometimes  inverted  the 
rating  between  two  materials. 

The  writer  found  that  Equation  (3)  held  fairly  well  for  metals  and 
simple  inorganic  compounds  as  diluents  of  TNT,  but  was  quite  inade¬ 
quate  for  dilution  with  air  or  a  plastic.3  The  available  data  on 
diluents  were  insufficient  to  test  the  Copp  6  Ubbeiohde  treatment. 
Instead,  an  expression  for  the  effect  of  entrainment  was  derived. 

This  was: 


Occurrence  of  a  phase  change  in  Si02  in  this  pressure  range  has 
been  verified.  See,  e.g.,  footnote  aa  of  Table  1. 
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D/Da  =  (1  +  C.5  yd)(l  -  0.58  yd) 

il+0.5yd  (1  -  3),  yd  small  (5) 

where  8  =  (pd  -  p^)/pe.  Equation  (5)  is  an  approximate  form  of 
Equation  (3)  showing  the  dependerce  on  yd.  When  this  equation  was 
modified  to  include  heating  the  diluent  by  compression,  it  became 

D/Da  =  1  +  yd(0.5  -.0.58  -  aot),  yd  small  (6) 

where  the  diluent  compressibility  ot  =  1  -  pd/pd’,  a  =  pj/(p-}  -  pe)*, 
pd '  is  the  density  of  the  shocked  diluent  and  Pj  is  the  C-J'density 
of  the  pure  -explosive.  The  approximate  linearity  of  D/Da  with  yd  for 
small  yd  was  substantiated  with  a  number  of  solid  diluents.3 
Figure  2  demonstrates  it  for  A1  and  NaCl.  . 

By  Equation  (5),  the  older  data3  indicate  a  compressibility  of 
12%  for  NaCl ;  the  data  of  Figure  2  (for  NaCl  coarser  than  50 Op),  16%. 
Interestingly  enough,  Copp  £  U'nbelohde  found  13  tc  15%  for  NaCl  of  an 
average  particle  size  of  about  ^2li.  When  so  many  variables  are 
involved,  the  value  indicated  for  a  does  net  seem  a  very  sensitive 
test  of  the  initial  assumptions. 

With  the  same  procedure  at  A  =  0.95,  the  NOL  data  (Figure  la) 

indicate  a(20p  Al)  =  0.073  whereas  from  Figure  3  data,  a(0.2p  Al)  = 

0.154.  In  other  words,  a(Al)  appears  sensitive  to  a  particle  size 
effect  whereas  «(NaCl)  does  not.**  The  simplest  explanation  is  that 
NaCl  acts  to  lower  D  chiefly  by  its  compressibility  arid,  because  of 
its  poor  thermal  conductivity,  is  not  much  affected  by  particle  size; 
whereas  metals,  which  are  good  thermal  conductors,  can  absorb  heat 
from  the  hot  detonation  product  gases.  For  this  mechanism  the  total 
amount  of  metal  surface  exposed  to  the  gases  and  hence  the  particle 

size  would  have  a  large  effect.  Thus  for  Al  and  other  metals 

Equation  (6),  which  omits  the  heat  absorption  mechanism,  is  probably 
incomplete  and  certainly  inadequate. 

If  cx  is  estimated  from  data  of  Figure  2b(A  =  C.80  instead  of  0.95) 
and  Equation  (6),  the  value  found  for  Ai  is  nearly  triple  that 
calculated  for  tne  less  porous  charge;  the  computed  a  value  for  NaCl, 
however,  changes  only  about  10%  with  the  same  change  in  porosity. 

This  again  suggests  the  need  to  supplement  Equation  (6)  with  some 
consideration  of  heat  absorption  by  the  diluent.  Such  an  effect 
will  be  small  at  any  porosity  for  nonconductors  suen  as  NaCl;  it 
evidently  increases  with  porosity  for  metals,  as  is  reasonable  in 
view  cf  the  increasing  accessibility  of  surfaces  for  thermal  transfer 
and  che  increase  in  detonation  reaction  time,  i,e.,  increased 
exposure  time  of  diluent. 

’•Derived  for  A  -  1,  a  =  k  +  1  where  k  is  the  adiabatic  exponent  for 
the  explosive  at  its  crystal  density.  Tno  value  of  a  used  for  TNT  in 
reference  (3)  was  4.22.  The  value  used  nere  is  3.85  in  agreement 
with  recent  Ruby  code  results  for  TNT. 

**See  also  reference  (3). 
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In  summary,  TNT/A1,  like  the  TNT/inert  mixtures  of  reference  (3), 
shows  a  linear  variation  of  D/Da  with  y^,  the  volume  fraction  of  the 
diluent.  There  is  evidence  from  both  particle  size  and  porosity 
effects  that  Equation  (6)  is  adequate  only  when  the  diluent  is  in  a 
coarse  form,  e.g.,  average  particle  size  greater  than  250p,  or  is  a 
thermal  insulator  or  both.  Otherwise,  it  requires  a  modification 
to  account  for  heat  transfer  to  the  diluent  particles;  this  transfer 
is  particularly  important  for  metals  (good  heat  conductors). 

Equation  (6)  serves  to  emphasize  the  necessity  of  making 
comparisons  at  equal  volume  instead  of  at  equal  weight  fraction  of 
the  diluents.  It  also  shows  how  two  diluents  of  the  same  density 
(same  f>  value)  can  have  different  effects  on  D  i>y  virtue  of  different 
compressibilities.  The  present  discussion  has  also  pointed  out  the 
possibility  of  different  thermal  properties  effecting  different 
amounts  of  lowering  of  the  detonation  velocity. 

Table  1  shows  what  are  considered  the  important  properties  of  the 
diluent  in  determining  its  effect  on  D.  It  includes  the  diluents 
Al,  NaCl,  W,  and  Si02 .  As  would  be  expected,  the  thermal  conductivity 
of  the  metals  is  two  orders  of  magnitude  greater  than  that  of  the 
non-metals.  Al  is  a  better  conductor  than  W  when  both  are  solids; 
with  liquid  Al  (m.p.  660°C),  the  relationship  is  reversed.  It  takes 
very  nearly  the  same  amount  of  heat  to  raise  equal  volumes  of  solid 
Al  and  W  to  the  same  temperature,  and  that  amount  is  about  30%  greater 
than  the  thermal  energy  required  to  heat  the  same  volume  of  NaCl  or 
Si02  to  the  same  temperature.  The  major  difference  in  energy  absorbed 
in  the  short  (detonation  reaction)  times  will  therefore  be  determined 
by  the  rate  at  which  heat  can  flow  in,  i.e.  ,  by  the  thermal  conduc¬ 
tivity  of  the  material.  Thus,  the  effect  of  equal  volumes  of  Al  and 
W  at  (nearly)  the  same  particle  size  distribution*  should  be  about 
the  same  as  compared  to  the  very  different  effects  of  diluents  such 
as  NaCl  and  SiG^, . 

No  attempt  wil)  be  made  here  to  develop  a  model  accounting  for 
the  effects  on  D  of  diluent  density,  compressibility,  thermal 
diffusivity,  and  charge  porosity.  Rather.  I  will  present  a  method 
of  treating  the  relatively  few  data  available  so  that  some  conclusions 
can  be  drawn  about  the  effect  of  different  diluents. 

To  compute  the  extent  of  the  lowering  of  D  from  the  data  measured 
over  a  range  in  A,  Equations  (1)  or  \2)  might  be  used  to  give 


D/Da 


a  2  -  f~ 
ae  *>e« 


r  •»  \ 


Where  the  subscript  e  designates  the  explosive  matrix  of  the  diluted 
mixture.  However,  the  form  of  Equation  (7)  is  a  bit  awkward:  it  is 
more  convenient  to  use  the  form 


*These  two  conditions  cannot  be  exactly  satisfied  simultaneously 
unless  the  diluents  have  equal  densities.  This  restriction  will  be 
considered  later. 
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D  =  AA 


n 


( 8  ) 


so  that 


D/Da  =  AAAn 


(9) 


where  A  =  A/Ae  -  1  and  An  =  (n  -  ne).  Since  A  is  the  value  of  D  at 
A  =  1,  the  approximation  S’  =  1  is  best  for  TNT/A1  mixtures  and 
poorest  for  TNT/NaCl,  70/30,  (See  Figures  1  and  3).  Thus,  An  is  a 
single  dominant  index  of  the  effect  of  the  diluent  at  ail  A. 


Figure  5  shows  that  Equation  (8 )  fits  the  NOL  D-A  data  for 
0  to  32%  A1  to  well  within  experimental  error  in  D.  In  this  work, 
the  parameters  A  and  a  were  evaluated  for  a  range  of  Po  =  1  to  2  g/cc. 
Both  Equation  (1)  and  Equation  (8)  are  valid  only  when  restricted  to 
a  specific  range  in  A.  The  range  chosen  for  evaluation  of  the 
parameters  is  not  necessarily  the  optimum;  hence  the  equations 
obtained  are  not  necessarily  the  best  fits.  Nevertheless ,  they  serve 
the  purpose  of  the  illustration  quite  adequately. 


Table  2  contains  the  parameters  of  Equations  (1)  and  (8)  for  the 
various  diluents.  In  Figure  6,  the  diluent  effect  on  D  is  approx¬ 
imated  by  An  as  a  function  of  yd* ,  the  diluent  volume  fraction  at 
voidless  density.  W  with  its  high  density  of  19  g/cc  has  a  very 
small  volume  fraction  at  15%  by  weight.  The  NOL  iata  for  TNT/2  Op  A1 
are  bracketed  by  reference  (6)  data  for  0.2u  and  270y  Al.  The  NaCl 
data  (reference  (5))  show,  despite  the  scatter  of  the  data,  the 
small  effect  of  that  diluent.  The  similarity  of  metal  diluents  is 
indicated  by  the  fact  that  a  reasonable  extrapolation  of  the  TNT/0.4p 
W  will  pass  very  near  the  single  point  for  TNT/0.2p  Al.  This 
possibility  can  be  explored  further  by  computing  comparable  values 
of  D/Da. 


Although  An_may  be  considered  the  dominant  parameter.  Equation  (9) 
which  includes  A,  must  be  used  to  compute  D/Da  if  inversions  in 
ratings  such  as  that  of  Figure  2  are  to  be  found.  Table  3  compares 
the  computed  results  with  the  ratio  values  of  the  experimental  results 
at  the  higher  porosity  where  the  change  in  D  is  larger.  It  can  be 
seen  that  for  A  =  0.78  -  0.85,  the  computed  values  agree  with  The 
experimental  to  better  than  2%.  (This  is  just  a  test  of  the  adequacy 
of  the  parameters  A,  n  of  Table  2  in  reproducing  the  experimental 
data  from  which  they  were  derived.)  The  D/Da  values  show  that  on  an 
equal  weight  basis  and  at  A  =  0.815  ±  0.035,  Al  decreases  D  more  than 
does  W.  However,  if  we  make  a  reasonable  interpolation  for  0.2jj  Al, 
as  shown  by  the  dashed  curve  of  Figure  6,  on  an  equal  volume  basis  at 
A  =  0.85,  W  decreases  D  more  than  does  Al .  To  be  sure,  the  difference 
is  not  of  significant  size,  but  the  calculated  D/Da  values  Lend  tu  be 
too  low  for  f.l  and  too  high  for  W.  So  the  difference  is  probably 
real  and  small. 


On  an  equal  weight  basis,  with  particle  sizes  of  0.2  and  0.4u, 
respectively,  for  Al  and  W.  the  number  of  Al  particles  is  56.5  times 
the  number  of  W  particles  and  the  total  surface  of  Al  is  14.1  times 
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that  of  the  W.  In  contrast,  on  an  equal  volume  basis  (say  yd*  =  0.015 
and  xai  =  0.0243,  x^  =  0.15),  the  number  of  A1  particles  is  only 
9.2  times  the  number  of  W  and  the  surface  area  of  A1  only  2.3  times 
that  of  W.  Insofar  as  heat  absorption  is  important,  there  is  good 
reason  for.  0.2y  A1  to  lower  D  more  than  0.4p  W  on  an  equal  weight 
basis. 

If  the  comparison  could  be  made  at  equal  y^*  and  equal  particle 
size,  the  number  of  A1  particles  and  total  A1  surface  would  still  be 
1.15  times  the  corresponding  W  values.  This  is  because  the  ratio  of 
the  number  of  particles  and  of  the  surface  area  of  two  diluents 
varies  inversely  as  the  ratio  of  the  voidless  densities  of  the  two 
mixtures.  (For  yd*  =  0.015  that  ratio  for  W  and  A1  in  TNT  is 
1.912/1.657  =  1.15.)  Hence,  only  when  the  diluents  have  equal 
density  is  it  possible  to  have  equal  yd*  and  an  equal  number  of  equal 
diameter  particles  at  the  same  time. 

Reference  (6)  did  report  on  silica  as  a  diluent  and  listed  its 
density  as  2.70  g/cc,  ir  other  words,  equal  that  of  A1 .  The  authors 
then  concluded  that  Al  must  react  endothermally  at  A  =  0.85  because 
it  lowered  D  more  than  did  the  inert  Si02 .  On  the  other  hand,  they 
considered  W  inert  because  of  its  high  m.p.  (3370°C)  and  with  that 
I  agree.  As  was  remarked  above,  Dremin  et  al.  were  aware  of  the 
necessity  for  comparisons  at  equal  A,  but  they  seemed  unaware  of  the 
necessity  for  comparisons  at  equal  volume  of  diluent.  On  this  basis, 
Al  seems  no  more  effective  in  lowering  D  than  the  inert  W  when  they 
are  in  comparable  physical  forms.  Moreover,  the  lesser  effect  of 
the  non-metals  can  be  explained  on  the  basis  of  their  low  thermal 
conductivity.  Comparison  of  the  results  on  an  equal  A,  equal  yd* 
and  comparable  6  basis  substantiates  the  conclusion  of  reference  (1) 
that  Al  does  not  react  in  time  :o  contribute  to  the  detonation  front 
of  aluminized  organic  explosives. 

Aluminized  RDX  and  TNETB 


The  binary  TNT  mixtures  were  discussed  at  some  length  in  the 
preceeding  section  because  some  systematic  data  were  available  for 
different  diluents.  This  is  not  the  case  for  other  organic  K.E. 

NOL  data  for  the  RDX/A1  and  for  the  TNETB/A1  series4  can  be  treated 
as  were  those  for  the  TNT/A1  series.  The  Al  used  was  atomized,  met 
Navy  specifications  and  had  been  sieved  through  a  No.  40  (420p) 
screen.  As  for  TNT,  the  velocities  were  measured  at  a  constant 
charge  diameter  of  5.08  cm.  Both  aluminized  series  showed  D-A 
patterns  comparable  to  those  of  TNT/A1  Figure  la,  i.e.,  increasing 
Al  decreases  D  and  increases  the  slope  b.  The  P-A  curves  fan  out 
from  A  ~  1.0;  hence  the  greater  the  porosity,  the  larger  the  decrease 
in  D  caused  by  a  given  amount  of  Al.  Because  of  its  small  size, 
the  decrease  is  therefore  hard  to  measure  near  voidless  density. 

The  da+"a  for  tnese  two  .  ories  show  more  scatter  than  those  for 
TNT/A1,  but  again  D/Da  seems  a  linear  function  of  Vq  (at  low  contents 
of  Al)  and  relative  effects  can  be  reversed  at  different  porosities, 
as  shown  in  Figure  2.  At  5%  porosity  (A  =  0.95),  a  given  volume  of 
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A1  decreases  D(TNETB)  more  than  D(RDX)  whereas  at  20%  porosity  the 
reverse  is  true. 


Table  4  contains  the  parameters  of  Equations  Cl)  and  (8)  for 
the  RDX  and  TNETB  series,  and  Figure  7  shows  the  plot  cf  An  vs  yd* 
for  all  three  series.  Two  points  from  reference  (6)  have  been 
included  to  show  that  the  effect  of  A1  particle  size  is  much  larger 
than  any  difference  introduced  by  varying  the  matrix  H.E.  Therefore, 
we  cannot  tell  how  much  the  location  of  the  curve  is  caused  by 
particle  size  effect  and  how  much  should  be  attributed  to  the  organic 
H.E.  It  is  interesting,  however,  that  a  plot  of  An  vs  the  ratio 
Al/0,  shown  in  Figure  8,  avoids  any  intersection  of  the  curves.  In 
Table  5 ,  the  approximate  lowering  of  the  detonation  velocity 
[D/Da  =  A An]  is  shown  at  two  porosities  for  each  of  the  three 
aluminized  explosives  at  the  same  values  of  Al/0.  At  these  two 
porosities  (10  and  50%),  the  amount  of  the  lowering  increases  in 
the  order:  TNT,  TNETB,  and  RDX  series. 

Effect  of  Aluminum  on  Detonation  Pressure 

Because  addition  of  A1  lowers  the  detonation  velocity,  it  must 
also  lower  the  detonation  pressure  (P).  From  the  usual  hydrodynamic 
relations. 


and  also 


P  =  PqD2 / (k+X)  =  APvD2 / (k+1) 


P  -  PoDu  =  APwDu 


(10) 


(11) 


where  PQ  is  charge  density,  k  is  the  adiaDatic  exponent,  and  u  is 
particle  velocity.  If  we  again  let  subscript  a  indicate  the  pure 
explosive  at  the  same  A  as  the  diluted,  then  from  Equation  (10) 


and  from  Equation  (11) 


P  =  £v  ke  +  i 
¥&  Pe  k+1  \2a/ 


p  uD 
v 


(12) 


Pa  peue^a 

If  the  effect  of  the  diluent  is  only  a  reduction  of  kinetic  energy  by 
entrainment,  equation  (3)  can  be  used  with  Equations  (12)  and  (13) 
to  give 


(P/Pa)en  =  ye-t (ke+l )/ (k+1) 


(P/Pa)ert 


■st-W'-iH 


(14) 


(15) 


1st 


_ _ _ ,_ _ _  _ _  „  _ ___  _ «.~**rT* 


NOLTR  72-62 


t£>  CD  CD  .3-  i_D 
CM  rH  LO  CO  CM 
COtDCDCOC"* 

O  O  O  O  rH  CM 


co  o  co  cn 
to  rH  cd  =* 
H  co  zt  zt 
*  •  •  • 
O  O  O  O  rH 


id  cm  zr  r-*  o 

r-  cn  o  lo  o 

O  CD  cn  O  CO  f-H 

r-i  CO  '  ^  t-i  O  r-\ 


CD  CNJ  CD  CO 

Is*  to  O)  o 

O  CD  CD  CD  CD 
rH  O  O  O  O 


CO 

c 

c 

O 

•H 

-P 

aJ 

a 

3 

0) 

cr 

<  03 

u 

3 

(D  tO  lO  CO  fO  J 
LO  CO  r-  rH  CD  CO 
CD  CD  CM  CD  CD  CO 

O  O  t — I  rH  rH  CO 


CO  lO  CD  O  CD  LO 
CO  CD  CD  lO  C^ 


co  co  co  co  r-  cd 


zi  r-  ^  cm  co 

CM  C-  00  CM  C'- 
N  O)  O  CM  rH 

O  O  rH  rH  CM 


CM  Zf  O  CD  CO 
tO  CO  CM  Z}-  zT 


o  CO  OO  CO  00 


.3- 

M 

C 

co 

z 

o 

3 

o 

M 

•rl 

N.  -» 

I — 1 

z: 

P 

E  x 

XI 

D 

m 

E  Q 

<0 

XJ 

3 

Eh 

< 

cr 

nJ 

w 

CD  CM  CM  O  CM  CO 
J-  JCNIOHCO 

cm  a-  o  co  to 

CD  00  O  CM  CO  zl" 
i — i  H  i — I  H 


tO  O  O  O  CM  00 
rHOOOCDr- 
co  rH  CD  rH  LO  O 

CM  O  H  co  zT  CD 

lilt 


ID  H  J-  J-  GO 
H  CD  O  O  J 

lO  CO  J-  in  H 

CD  00  CD  CM 


o  o  o  co 

ZT  CM  O  O  CO 
CD  co  CM  CM  C) 
•  •  •  • 
HOOHJ“ 


CD 

OO  CO  CM  00  O 
cd  zr  cm  o  o 
O  rH  CM  CO  zr 

o  o  o  o  o  o 


CM  J- 

CO  rH  O  OO 
CD  ZT  CM  o> 
O  rH  CM  CO 

o  o  o  o  o 


o  o  o  o  o 

rH  CM  CO  .-3-  LO 


o  o  o  o 

H  CM  CO  LO 


o  o  o  o  o  o 


o  o  o  o  o 


(0 

to  >*  a 
o)  -h  a 

rH  -H  \ 
T3  U)  be 
•H  C  ^ 
O  Q)  > 
>  Q  a 


CM 

O  CD  CO  O  00  CD 
00  CO  CD  O  O  H 

rH  rH  rH  CM  CM  CM 


00  zT  H  00  tO 
C''*  CO  CD  CD  rH 


rH  rH  'H  rH  CM 


CO  CO  CO  tO  rH 
tOzTOOrHr^ 
rH  CO  tO  CD  CO 

O  CO  CO  CO  CO  r-i 


.1  (£>  tO  NJ 
H  lO  CO  CM 
rH  CM  zJ-  O 
•  •  •  • 
O  O  O  O  rH 


CHANGE  IN  EXPONENT,  An 


>-S-^”'*>vy-J 


~T.-:-K^ 


NOLTR  72-62 


FIG.  7  CHANGE  IN  EXPONENT  An  EFFECTED  BY  DIFFERENT 
MATRICES  IN  HE/A!  MIXTURES 
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Table  5 
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Dremin  et  al .6  measured  u^j  and  D  of  TNT  and  85/15,  TNT/diluent 
at  96.5  and  85%  TMD.  They  found  that  much  of  their  data  was  well 
expressed  by 

u  =  ueAp_/£pv 

6  (16) 

-  Ug  Pg /  Py 

where  u  and  ue  are,  respectively,  the  C-J  particle  velocities  of  the 
mixture  and  of  the  pure  explosive.  This  expression  in  Equation  (13) 
results  in 

(P/Pa)D  =  D/Da,  (17) 

an  unexpected  result  because  it  indicates  that  the  aluminum  lowers 
the  pressure  by  the  same  percentage  that  it  lowers  D  rather  than 
D2.  The  few  pairs  available  (HBX-1  £  KBX-3,  and  Comp  B  and  H-6)  do 
not  validate  Equation  (17),  although  those  data  of  reference  (6)  for 
which  Equation  (16)  holds  do  conform  to  Equation  (17). 

Equation  (16)  holds  for  the  aluminized  charges  (0.2  and  270y  Al) 
at  96.5%  TMD  and  to  within  6%  for  the  charges  (80  and  270p  Al)  at 
80%  TMD.  It  is  off  by  21%  for  0.2y  Al  at  the  lower  %  TMD.  The 
authors  attribute  the  6%  discrepancy  to  heat  absorption;  they 
computed  a  decrease  of  up  to  5%  for  this  mechanism.  For  the  very 
fine  Al  at  the  higher  porosity  they  believe  the  large  decrease  is 
caused  by  endotnermal  reaction  between  the  Al  and  explosive .  From 
qualitative  arguments  above,  heat  loss  by  absorption,  compression, 
and  entrainment  seem  sufficient  to  explain  the  observed  effects; 
semi-quantitative  considerations  are  given  in  the  Appendix.  As  I 
pointed  out  earlier,  the  possibility  of  endothermic  reaction  is  still 
unresolved. 

The  fact  that  Equation  (16)  holds  for  aluminized  TNT  at  high  TMD 
suggests  the  possibility  of  better  approximating  pressure  from 
measured  D  values  in  cast  aluminized  explosives.  The  usual  approx¬ 
imation  P  =  0.25poD2  generally  gives  too  high  a  value  for  aluminized 
H.E.  because  k  is  greater  than  3  and  (k+1)-1  <  0.25.  On  the  other 
hand,  the  usual  approximation  fits  the  reference  (6)  data  for  high  A, 
aluminized  TNT  to  1.5%  or  better  than  it  fits  the  data  for  pure  TNT 
(to  5.8%).  The  approximations 


P  =  0 . 25ApeD2 
P  =  0 . 2  5 APPWD2 


(18a) 

(18b): 


as  well  as  the  conventional  one  have  been  used  for  the  comparisons  of 
Table  6.  In  view  of  the  uncertainty  in  the  measured  values  of  P 
and  also  of  the  D  value  for  HBX-3,  the  results  are  inconclusive. 
Equation  (18b)  seems  to  be  a  slightly  better  approximation  for  the 
four  explosives  of  Table  6  than  Equation  (18a);  both  seem  better 

*Pew  is  voidless  density  of  the  matrix  including  the  wax. 
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xhan  the  usual  PQD2/4.  If  this  proves  to  be  the  case  when  better 
data  are  obtained,  it  might  still  be  restricted  to  this  particular 
type  of  explosive  (aluminized  Comp  B)  within  a  restricted  range  of 
aluminum  content. 

Lyubov  et  al.12  nave  derived  an  expression  for  the  detonation 
pressure  of  explosives  containing  an  inert  additive;  they  consider 
A',  inert.  In  our  notation,  their  relationship  is 


or  (19) 

P  =  P0D2ye/ (k+1 ) 

This  is  the  usual  hydrodynamic  relation  [Equation  (10)]  with  a 
correction  factor.  The  authors  claim  that  detonation  pressures 
calculated  according  to  Equation  (19)  for  mixtures  of  PETN  with 
Al  and  Mg  as  well  as  of  TNT  with  Al,  Si02 ,  and  NaCl  (data  of 
reference  (6))  containing  up  to  50%  inert  additive  agree  with  the 
experimental  values  within  10%.  Their  tabulation  shows  that  they 
used  a  constant  value  of  k  =  2.80  for  the  calculations  regardless 
of  matrix,  diluent  content,  and  charge  porosity.  Hence,  Equation  (19) 
can  hardly  be  considered  an  established  relationship. 

In  summary,  addition  of  Al  to  organic  explosives  lowers  the 
detonation  pressure.  There  is  no  convenient  approximation  available 
with  which  to  estimate  the  extent  of  the  decrease. 

SUMMARY 

A  review  and  examination  of  available  data  for  aluminized  organic 
explosives  has  shown  that : 

1.  In  all  probability,  the  aluminum  does  not  react  in  the 
detonation  front  but  sometime  later. 

2.  Al  as  a  ciluent  lowers  D.  The  lowering  is  a  linear  function 
of  the  volume  fraction  of  Al  present;  it  also  depends  on  the  dens5.ty 
and  compressibility  of  Al  as  well  as  or.  its  particle  size  and  shape  — 
probably  on  the  total  surface  and  heat  conductivity. 

3.  The  lowering  of  D  by  dilution  with  Al  also  depends  on  the 
particular  H.E.  used.  The  same  amount  and  form  of  Al  has  a  larger 
effect  (expressed  as  %D)  on  explosives  of  initially  higher  D  values. 
Thus  the  effect  decreases  in  the  order  RDX,  TNETB,  TNT. 

4.  Aluminum  as  a  diluent  lowers  P. 

5.  The  decrease  in  D  and  P  is  most  obvious  at  high  porosity 
(low  A)  and  high  Al  content;  the  former  is  the  more  important  factor. 
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SYMBOLS  AND  DEFINITIONS  USED 

Component  i  of  mixture  (also  ideal  value  of  D) 

Explosive  component 
Solid  diluent 
Air 

Voidless 
Initial  value 

Reference  value  D  or  P  of  pure  explosive  at  corresponding 
porosity. 

Weight  fraction 
Volume  fraction 
Specific  volume 
Bulk  density  of  charge 


A  =  PQ/pv  =  i  -  ya 


Pv 


(Solids  \-1 

?  *i'Pil 


y±  =  Apvxi/pi  =  p0X;i/pi 

a  =  Compressibility  of  diluent  =  1  -  P^/P^’  =  Cv^-v^' )/v^  where 
prime  indicates  compression  in  detonation  front. 

6  -  ( Pc[/  Pe )  ~  1 

a  =  1/ (1  -  pe/pj)  =  k  +  1 

k  Adiabatic  exponent 

a(u),  b(u)  Parameters  of  D  =  a  +  tA 

A,n  Parameters  of  D  =  AAn 
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Appendix  A 


ADDITIONAL 


CONSIDERATIONS  OF  ROLE  OF  ALUMINUM  IN  H.E./A1  MIXTURES 


Reference  (6)  does  not  give  the  calculations  made  for  heat 
absorption  by  A1  and  the  consequent  effect  on  D.  In  view  of  our 
disagreement  with  the  conclusions  of  reference  (6),  it  seems  worth¬ 
while  to  present  such  estimates  as  we  can  carry  out  with  our  present 
knowledge  of  the  characteristics  of  Al. 

Lees13  has  studied  the  melting  point  of  Al  under  pressure  over  the 
range  0  to  40  kbar.  His  data  show  the  linear  relation 


Tm  =  660  +  6.58P  (Al) 

where  the  units  of  T,  and  P  are,  respectively,  °C  and  kbar. 

Krieger11*  has  studied  the  thermodynamics  of  the  aiuminum/aluminum 
vapor  system,  and  constructed  the  vapor  pressure  curve.  From  that 
curve,  the  boiling  point  of  Al  is  4300°K  at  0.1  kbar,  and  greater 
than  600G°K  at  1  kbar  where  it  is  still  rising  rapidly.  Hence  the 
Al  may  melt  in  the  detonation  zone,  but  it  will  not  vaporize. 
Consequently,  any  diffusion  reaction  with  the  H.E.  detonation 
products  will  be  restricted  to  the  surface  of  the  Al  melt,  much  as 
it  is  in  the  case  of  aluminized  rocket  propellants. 


To  consider  the  effect  of  heat  absorption  let  us  assume  that 
the  aluminum  is  exposed  to  the  C-J  temperature  (Tj )  at  the  C-J 
pressure  Pj  for  an  interval  equal  to  the  detonation  reaction  time  x. 

The  detonation  reaction  zone,  in  the  one  dimensional  case,  is 
bounded  by  the  leading  von  Neumann  (vN)  shock  plane  and  the 
following  C-J  plane  where  the  reaction  for  C-H-N-0  explosives  is 
essentially  complete.  From  the  vN  plane  to  the  C-J  plane,  the 
pressure  decreases  and  the  temperature  increases.  Hence  Tj  is 
higher  and  Pj  is  lower  than  the  average  values  over  the  interval  x. 
However,  there  is  no  general  agreement  on  the  correct  value  of  Pj 
(measured  or  computed)  and  even  less  on  the  value  of  Tj ,  which  may 
vary  by  several  thousand  degrees  according  to  the  equation  of  state 
used  for  the  gas  products.  Under  these  circumstances,  there  seems 
little  point  in  trying  to  pick  average  values;  Table  Al  contains  the 
Pj ,  Tj  values  for  TNT  (and  their  sources)  chosen  to  estimate  the 
order  of  magnitude  of  the  heat  effecx.  Also  tabulated  are  measured 
values  of  the  detonation  reaction  time  x.  (This  is  not  to  be 
confused  with  a  total  reaction  time;  it  merely  indicates  the  interval 
beyond  which  reaction  energy  no  longer  contributes  to  propagation 
of  the  detonation  front.)  The  values  chosen  for  the  measured  Pj 
and  the  computed  Tj  are  the  lower  of  those  available  in  the  literature; 
the  value  of  x  is  an  intermediate  of  those  measured. 
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In  addition  to  the  above  values,  we  need  information  about  Al. 

Its  melting  point  at  the  higher  pressures  is  estimated  by  extra¬ 
polation  of  Equation  (Al).  The  thermal  conductivity,  specific  heat, 
and  density  of  Al  will  all  vary  with  temperature  and  pressure  and 
so  to  will  its  thermal  diffusivity.  Because  data  are  lacking,  the 
value  k  =  0.94  cm2  sec-1  computed  from  the  values  of  k/spv  of 
Table  1  is  used  for  diffusivity.  The  necessary  parameters  Kt/a2 , 
where  a  is  the  average  particle  radius,  are  then  those  tabulated. 

Solution  of  the  heat  conduction  equation  for  spheres  with  the 
surface  temperature  Tj  for  the  time  t,  shows  that  the  entire  particle 
of  the  0.2y  Al  is  raised  to  a  temperature  above  its  melting  point 
(in  fact,  it  would  be  raised  to  Tj  in  time  t  if  its  surface  could 
be  maintained  at  T-j )  whereas  only  the  outer  layers  of  the  larger 
(80  and  27  0y)  particles  can  be  heated  above  the  melting  point. 

It  follows  that  with  sufficient  time  the  0.2y  particles  could  melt 
completely  whereas  the  larger  particles  would  show  only  a  thin  layer 
of  melt  on  the  surface.  These  data  illustrate  the  particle  size 
effect  on  heat  absorption  rather  more  strongly  than  the  A  effect, 
but  that  too  is  evident. 

Although  there  have  undoubtedly  been  many  thermodynamic- 
hydrodynamic  calculations  of  the  detonation  of  aluminized  organic 
explosives,  very  few  have  been  published  in  accessible  sources. 

Mader16  used  the  LASL  code  for  one  aluminized  H.E.,  the  composite, 

RDX/TNT/Al/wax,  44/32.2/19.8/4  at  A  =  0.9945.  Not  only  is  this  at 
a  compaction  where  little  or  no  effect  of  Al  on  D  is  expected,  but 
he  failed  to  run  the  computation  on  the  non-aluminized  matrix,  j 

RDX/TNT/wax,  54.86/40.15/4.99.  With  the  assumption  of  complete  I 

reaction  of  Al  to  form  Al?03(s),  he  computed  D  about  3%  below  the  f 

experimental  value  of  7.53  mm/ y sec  and  P,  3%  above  the  experimental  l 

value  of  230  kbar.  The  Kamlet  empirical  method18,  which  approximates  1 

very  well  the  current  RUBY  code  values  for  simple  C-Ii-N-0  explosives,  ] 

can  be  used  to  estimate  the  parameters  for  the  non-aluminized  matrix.  = 

This  was  done  and  the  wax  was  treated  as  a  reacting  component  to  • 

obtain  the  values  of  7.60  mm/ysec  and  245  kbar  respectively  for  D  < 

and  P  of  the  matrix.  These  values  are,  within  the  error  of  the  $ 

method18,  the  same  as  those  measured  on  the  aluminized  charge  at  j 

A  =  0.9945.  Hence  they  confirm  the  insignificant  effect  of  Al  at  i 

high  compaction  (relative  to  its  matrix  at  the  same  A)  but  not  the  ] 

reason  for  this  effect.  Thus  at  A  ~  1  Mader  computes  the  same  ; 

result  assuming  reaction  of  the  Al  that  Equation  (6)  would  give  for  ] 

chemically  inert  Al  that  is  compressed  and  entrained.  f 

One  other  set  of  computations  has  been  carried  out  by  Chaiken19  | 

on  the  RUBY  code  for  an  RDX  series  at  A  =  1,  again  where  the  effect  [ 

of  Al  should  be  insignificant.  His  results  are  shown  in  Table  A-2 ,  j 

but  they  are  very  puzzling.  The  present  LASL  code  is  very  similar  .  j 

to  the  RUBY  code.  There  seems  to  be  no  reason  for  Chaikin’s  results  i 

for  pure  RDX  to  be  so  low;  the  LASL  values  agree  with  the  experimental  j 

ones.  If  we  assume  an  error  in  input  for  that  explosive  only,  we  can  | 

compare  the  RDX/Ai  results  to  the  LASL  RDX  results.  That  comparison  '  | 

says  that  the  Al  has  practically  no  effect  on  D  and  P  at  A  =  1  when  \ 

Al  is  assumed  to  react.  The  effect  of  melting  the  Al  but  not  1 
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Table  A1 

TNT  AND  TNT/A1,  85/15 


TNT 


Reference 

A=0 . 9  64 

Pq=1 . 59 
g/cc 

A=0 . 85 
Po=1.40 
g/cc 

D(mm/ysec) 

15 

6.94 

6.34 

Pj (kbar) 

15 

175 

130 

Tj(°K) 

16 

2978 

3100 

ifysec) 

17,15 

0.14 

0.29 

Al,Tm(°C) 

E.i.  (Al) 

1811 

1515 

Tm  (  °  C  )  / 

Tj (°C) 

0.67 

0.54 

TNT/A1,  85/15 

a(cm) 

icx/a2  at  T>Tm* 

KT/a2 

at  T>Tm* 

0.01350 

0.00073  ~ly 
outer 
layer 

0.0015 

2 . 7y 

OUt€  ■ 

layer 

0.00400 

0.0082  2 . 5y 

outer 
layer 

0.017 

5 . 2y 
layer 

0.00001 

1330  whole 

2670 

whole 

sphere  sphere 


*From  Cars low  S  Jaeger,  "Conduction  of  Heat  in  Solids"  Oxford  Press, 
Oxford  (1959)  pp  234  S  102 
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allowing  it  to  react  seems  to  be  a  4%  increase  in  D.  This  is  hard 
to  believe  unless  4%  is  typical  of  the  error  in  the  program.  In 
view  of  the  reported  result  for  RDX,  it  seems  that  results  for  the 
whole  series  must  be  considered  unacceptable. 


In  earlier  work20  we  ha/e  shown  that  Kamlet's  method  does  not 
work  well  for  aluminized  explosives  (A1  assumed  reacting),  partic¬ 
ularly  as  A  decreases.  If,  however,  we  assume  the  A1  inert  and 
the  relationship 


D  =  n^Q^U.Ol  +  1.31pvA)* 


(A2  ) 


correct,  we  can  use  it  to  obtain  an  estimate  of  how  much  effect  heat 
absorption  could  have  at  A  =  0.85.  When  the  A1  is  0.2p  diameter, 
data  of  Table  1  indicate  that  the  particles  can  be  melted  and  the 
melt  heated  to  Tj  in  the  interval  r.  (This  is  an  upper  limit  to  the 
effect  because  the  product  gas  environment  could  not  maintain  a 
temperature  as  high  as  Tj  at  the  surface  while  it  is  losing  heat 
to  the  Al.)  Hence  to  apply  Equation  (A2)  to  TNT/A1 ,  85/15,  the 
pure  TNT  values  of  n  and  Q  must  be  reduced  by  the  factor  0.85;  M 
stays  the  same.  In  addition,  Q  must  be  further  reduced  by  heating 
the  Al  to  its  melting  point,  melting  it,  and  heating  the  melt  to 
Tj.  Thus  for  TNT,  Q  =  1282;  for  TNT/A1 ,  85/15,  Q  =  976.  In 
Equation  (A2)  this  gives  D  =  5.63  mm/ usee  as  compared  to  the  value 
of  5.75  mm/psec  measured  at  A  =  0.85.  Thus  heat  absorption  alone 
could  account  for  the  lowering  of  D  observed  at  this  compaction. 

The  example  is  very  probably  an  over-estimate  of  this  factor  as 
entrainment  and  compression  are  also  contributing. 


*n  =  moles  of  gas/gH.E.,  M 
Q  =  chemical  energy/gH.E. 


=  average  molecule  weight  of  gas  and 
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